CONCENTRATION OF THE RATIO BETWEEN THE GEOMETRIC AND 

ARITHMETIC MEANS 

do . 

O . J- M. ALDAZ 

o 

Abstract. We explore the concentration properties of the ratio between the geometric mean 
and the arithmetic mean, showing that for certain sequences of weights one does obtain 
concentration, around a value that depends on the sequence. 

o 
m 



1. Introduction. 



This paper is motivated by the article [GluMi] by E. Gluskin and V. Milman, who consid- 

r^ ! ered the variant YYi=i IVil^^"' — a/?^^^ XliLi vf °^ ^^^ AM-GM inequahty in the equal weights 

c^ ! case. Roughly speaking, they showed that the ratio YYi=i \yi^^^ I \/'^^^ XliLi v1 i^ bounded 

below by 0.394 asymptotically in n and with high probability, where probability refers to 

Haar measure on the euclidean unit sphere §2"^- Thus, the geometric and arithmetic means 

are comparable quantities on "large sets", provided n is large. 

^ ! For the most part, we focus on the standard AM-GM inequality HiLi klT " — X^iLi "^i.nkk 

^^ \ where aj,„ > and X]r=i '^«." ~ -'-• ^^ shall see that in some cases the concentration of 

00 ■ measure phenomenon does take place: For certain sequences of weights, which include the 

. . equal weights case, the GM-AM ratio is "almost constant", with probability approaching 1 

L^ ! as n — > 00. Here probability means normalized area over a suitable unit sphere (a level set of 

(X) \ the arithmetic mean /(x) = ^YTi=\ '^j.nkli) which depends on the weights. We prove that the 

O ■ GM-AM ratio cannot concentrate around values larger than e""^ ~ 0.5615, where 7 is Euler's 

constant. On the other hand, for every t G [0, e~'^\ it is possible to find a family of weights 

k> \ such that concentration takes place around t. In particular, when all weights are equal (to 

?H \ 1/n for ra = 2, 3, . . . ) there is concentration around e~"' , a result to some extent anticipated 

in |Glt Theorem 5.1]. 

For completeness, we also study Gluskin and Milman's equal weights modification of the 



GM-AM ratio: YYi=i Iz/il "/y ^^^ X^ILi vf- Here the natural choice of probability is the uni- 
form measure P^""*^ on §2 ""^ := {||y||2 = 1}, where ||?/||2 := \/Yll=i vf i^ the P^ norm on W^. 
For n » 1 and on Sg"^) the preceding ratio concentrates around -\/2 exp [r'(l/2)/(2r(l/2))] '^ 
0.5298. 

The method of proof used here is essentially the same as in |GluMi] : Compute the s- 
moment of HiLi I^IT" ^^*^ then use Chebyshev's inequality. The differences lie in the fact 
that we consider the usual AM-GM inequality rather than a variant of it, we do not restrict 
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ourselves to the equal weights case, and we optimize over s to obtain concentration results 
rather than lower bounds. 

The statistical properties of the ratio between the geometric and the arithmetic means, 
in addition to its mathematical interest, may be relevant in those areas where this ratio is 
used as a measure of homogeneity or information, such as, for instance, radar imaging (cf. 
IBMTEj . rWo] l 

2. Definitions and results. 

Given n >2, weights ai^n > with J2i=i ^i,n = 1) ^^'i ^ = {^ly ■ ■ ■ i ^n) G I^", the AM-GM 
inequality states that 



(2.0.1) nN"'"^E 



^2,n 1*^ I 
1 



Assume x ^ 0. Since the middle term of 






(2.0.2) < i[r '', < 1 

is a positive homogeneous function of degree zero, it is constant on the rays tv, where t G M 
and V G M". Thus, we may fix any suitable t and then select v G {}21=i^i,n\^i\ — "^^ 
random", that is, uniformly over this level set. Observe that the notion of random choice we 
use depends on the sequence a„ := (ai,n, • • • , «n,n)- Let us renormalize the weights by setting 

(2.0.3) ai^n ■= nai^n- 

Define the unit ball in R" associated to a„ := (ai,„, . . . , a„^„) by B"^ := {J2^=i (^i,n\xi\ < 1}, 
and denote the corresponding norm and unit sphere by || ■ ||a^ and S"^^ := {Y^^=i 0'i,n\xi\ = 1} 
respectively. In the case of equal weights, aj^„ = 1, || • ||a„ is just the £i-norm || ■ ||i, B"^ the 
cross-poly tope B", and S"^"*^ the £i-unit sphere S"~^. For arbitrary a„, B"^ and §"^^ are 
simply linear images of B", and S"~^. 

To select v G (X^iLi '^«.'il^«l ~ ^^ random, by homogeneity we can just set t = 1, and 
then use the normalization of area on S"^^. Denote by 7i^ the s-dimensional Hausdorff 
measure on M". We adopt the convention that "W is defined via the Euclidean i^ metric, 
and normalized by the factor tx^I"^ IY{\ + s/2), so if s = rf is a positive natural number, the 
cube of sidelength 1 has Hausdorff d measure 1. A consistent choice of normalizing factor is 
required, for instance, to use Fubini's Theorem, or more generally, the coarea formula. We 
denote by -P^~^ the uniform probability on S"~^, that is -P^~^ = c?i"~^, where c is chosen so 
that P"^~^(S"~^) = 1. Occasionally, and for simplicity, we will utilize absolute value signs to 
denote volumes and areas. 

The first three results deal with families of uniformly bounded renormalized weights. The 
fourth considers a sequence with some weights approaching cxd, which yields concentration 
around zero. 
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Theorem 2.1. Suppose that there exist an M such that ai^n < M for all n > 2 and all 
i = 1, . . . ,n. Let 7 be Euler's constant, let k,e > 0, and let P"~^ be the uniform probability 
on ^an^- Then there exists an N = N{k,e) such that for every n > N, the ratio of the 
geometric mean over the arithmetic mean on S^^^, given by T^HiLi I^jI"''"? satisfies 

(2-1-1) p::-' 1(1 - ^)^ < ^ n \^r- < (i + e)eA >^-^,- 

In the equal weights case we can set M = 1, obtaining concentration around e~^ ~ 0.5615. 

Corollary 2.2. With the same notation as in the preceding theorem, suppose additionally 
that ftjn = 1/n for i = 1, . . . ,n. Then there exists an N = N{k,e) such that for every 
n > N, on S"^"*^ we have 

(2.2.1) Pr' I (1 - s)e'^ < nH |x,|i/" < (1 + e)e-^ \>l-l.. 

Remark 2.3. The bounds given in Theorem 12.11 are optimal in the sense that the upper 
bound for concentration is achieved in the equal weights case, and it is also possible to obtain 
concentration around values arbitrarily close to e~^ /M for every M > 1, cf. Remark 12.61 
below. 

Remark 2.4. When studying the concentration of the GM-AM ratios, in addition to the 
uniform probabilities on the spheres S""^, there are other reasonable notions of random 
choice of a vector. Suppose, for instance, that we are in the equal weights case ai^„ = 1, with 
§^~i = §"^^. Denote this particular GM-AM ratio by ri„. We might, for example, select 
points from the whole space M", according the exponential density 2~"e~"^"^ on M", or the 
uniform density over B"(r). Of course, these probabilities give more weight to small vectors 
than to large ones. But since the GM-AM ratio is homogeneous of degree zero, this makes 
no difference: The averages of the GM-AM ratio are equal over all spheres §"~^(r) centered 
at zero with radius r > 0, so by the coarea formula, the same answer is obtained when using 
the exponential probability over M" or the uniform probability over B"(r). 

After having proved Corollary 12. 2[ I came across a related, large sample statistics result in 
the literature, cf. [GH Theorem 5.1], which in the special case k = 1 states that ri^„ — > e~'^ 
with probability 1 as n ^ 00, where the independent, identically distributed random variables 
Xi take values in [0, 00) according to an exponential distribution. Thus, probability refers to 
the corresponding product probability in the countably infinite product of positive semiaxes. 
Informally, this result has the same content as Corollary 12.21 If one chooses a large quantity 
of non-negative numbers at random, then with very high probability the GM-AM ratio will 
be very close to e~^. Formally however. Corollary 12.21 seems to be stronger, since by rewriting 
it in terms of the exponential distribution, it yields the above result, while the other direction 
does not follow from the statement of [Gh Theorem 5.1] (a problem is that when defining 
a probability on an infinite product many new measure zero sets may appear, which might 
have large outer measure from the n dimensional viewpoint). The proof of |Gll Theorem 5.1] 
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is not provided there, so I do not know whether it might yield something similar to Corollary 

Concentration may occur around levels different from e~^ . In fact, given any t G (0, e~'^] it 
is possible to find a uniformly bounded (by a fixed M = M{t)) sequence of weights such that 
the ratio of the geometric mean over the arithmetic mean concentrates around t. Observe in 

M-l 

the next result that M J^^+i maps {1 < M) onto (0, 1]. 

Theorem 2.5. Let 1 < M and let L » M . For n > L define the sequence {aj^„ : i = 
1, . . . ,n} of weights as follows: aj^n = M whenever j < n/{M + 1), aj^n = l/M whenever 
j > 1 + n/{M + 1), and if j is the least integer strictly larger than n/{M + 1), then aj^n £ 
[1/M, M] is chosen so Yl^=i^i,n = ^- Then there exists an N = N{k,e) such that for every 
n> N, the GM-AM ratio on E)"^'^ satisfies 

(1 - e)A^ < nUlx.r- < (1 + e) Attt > 1 - -• 

Remark 2.6. Note that 1/M m+i is close to 1/M provided M is large. To show optimality of 
the lower bound in Theorem 12. II for small values of M > 1, we can modify the finite sequences 
given in the preceding theorem by using 1/M^ as the small weights (with j » 1) instead of 
1/M, and a few more large weights M, in order to add up to n. Then the argument proceeds 
as in the proof of Theorem 12.51 

Suppose next that the largest weights for the ri-th averages are given by a function f{n) < n 
with lim„ f{n) = oo. Arguing as in Theorem [23] we obtain concentration of the GM-AM ratio 
around 0. 

Theorem 2.7. Let f{n) < n satisfy lim„/(?7,) = oo, and define {aj^„ : i = l,...,n} by 
0-3,11 = f{n) whenever j < n/{f{n) + 1), Oj^n = ^/f{n) whenever j > 1 + n/{f{n) + 1), and 
if i is the least integer strictly larger than n/{f{n) + 1), then aj^n ^ [l//(^)) /('^)] ^^ chosen 
so J2^=i ^i,n = ^- Then there exists an N = N{k, e) such that for every n> N , the GM-AM 
ratio on S"~"^ satisfies 

(2.7.1) pvUii 



xJ"'-" < e 




For completeness, we consider the equal weights variant of the AM-GM inequality 

n 

(2.7.2) ni!''l""< 



\ 



n 

n ^-^ 

1=1 



studied in |GluMi] . The arguments are essentially the same as in the previous cases, save 
that now the appropriate notion of random selection is given by the uniform probability P^~^ 
on the euclidean unit sphere §2^^ = {\\y\\2 = 1} (alternatively one might use, for instance. 
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the standard gaussian measure on M" to obtain the same conclusion, by the reasons given in 
remarkEaD. On §r\ Uli M"V \[lTJUyf = V^K=i bd^^"- 

Theorem 2.8. Let k,e > and let P^^^ be the uniform probability on E>2^'^ . Then there 

exists an N = N{k, e) such that for every n > N, 

(2.8.1) 

- ^)^-P &j] < Vr^fl |y.|'/" < (1 + .)V2exp (^] } > 1 - i, 




W 



Remark 2.9. V^exp (^fM) ~ 0.5298. 



Remark 2.10. Since \ ■^Yll=iy'i — ^ X]r=i ?/« ^^ Jensen's inequahty, Theorem 12.81 shows 
that with arbitrarily high probability 

(2.10.1) r^,n{x) := V^^n 1 , > (1 " ^) V2exp ' ^'^ 



^ELtk.l -^ ^'-^^^'^{2T{1) 

on the euclidean unit sphere §2"^; provided n is sufficiently large. Analogous remarks can 
be made about the behavior of the GM-AM ratio obtained from (12.7.21) on §"~^, by using 
Corollary O 

Remark 2.11. Often concentration is a consequence of some type of uniform Lipschitz 
behavior (cf., for instance, Levy concentration theorem in 3|.19 p. 142 of [GroJ ) . This is not 
the case here. Consider the equal weights ratio ri „, for instance. If we take as our sample 
space either M" or B"(r), then ri „ is not even continuous, regardless of how rin{0) is defined. 
On §"~^ the Lipschitz constant of ri„ depends on n: Set x = (1/n, 1/n, 1/n. . . , 1/n) and 
y = (0, 2/n, 1/n, . . . , l/n). Then ri^„(a;) — ri^niy) = 1 — = 1, while ||a; — y\\i = 2/n. 

3. Lemmas and proofs. 

Let us recall the coarea formula (for additional information we refer the reader to [E], PP- 
248-250, or [EG], pp. 117-119): 



(3.0.1) / g{x)\Jf{x)\dx= / / g{x)dW^{x)dt. 

Here / is assumed to be Lipschitz, and \Jf{x)\ := -^/det df{x)df{xy denotes the modulus 
of the Jacobian. We also remind the reader about some well known facts regarding the F 
function (cf. |Lu] . for instance): i) It admits the asymptotic expansion 

(3.0.2) T{z) = e-'z'-'/^V2^ (l + — + 0{z^^] 

ii) r(l/2) = vr^/^. iii) r'(l) = —7 (Euler's constant), iv) The logarithmic derivative of the F 
function is an analytic function on C \ {0, —1, —2, . . . }; in particular, it is continuous there. 
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The following notation, introduced in the preceding section, is recalled here: aj,„ > 0, 
^"=1 ttj.n = 1, and ai^n = nai^n- Without loss of generality we assume that the weights are 
arranged in decreasing order: ai^„ > a2,n > ■ ■ ■ > cin,n- The associated unit sphere and the 
uniform probability on it are denoted by S"~^ and Pa~^ respectively, and the norm, by 



an- 
n-1 



Lemma 3.1. Let s eM. satisfy l + sai„ > 0. Then the expectation ofYYi^i {xil""'-"^ overE)]^ 

is given by 

(3.1.1) 

E (f] \xT"A ■■= I if] i-.r-'-l dpr\x) = -p^, fn ^-^^ 

By renormalization of the weights, ai_„ > 1. The restriction s > — l/ai.„ is of no con- 
sequence to us since we will be using this result for s < (to obtain lower bounds) when 
ai,n < M and s ^0 (but s » l/n). 

Proof. It is well known and easy to check that the volume of the (M", || ■ ||i)-unit ball is \M^\ = 
2^/n\. Let T„ be the linear transformation satisfying T„(S"~^) = S"~^ = {X]r=i '^i,nki| = !}• 
Then detT„ = nr=i ^i~n- Now for every a; G M" \ U^^i{xi ^ 0}, the function f(x) = \\x\\a„ is 
differentiable, and df {x)df {x)* is a 1 x 1 matrix, so 



\Jf{x)\ = ^det{df{x)df{xy) 



E2 _ II II 



a.e. on M". Set g{x) = 1/| J/(x)| in (13.0.11) . and denote by §"^ ^{p) the sphere centered at 
of radius p (when p = 1 we usually omit it). Then 

^" |T„(B^)| = |B" 1= / dx= [' [ J—dn--\x)dp 



1 T-rn l^"V"'i/l I an I I I I II II 

^'llj=l«i,ri JBJJ„ Jo Js2nHp) Il'^™ll2 



Sn— 1| /•! |<Sn— 1| 

Pnib Jo ^llflnlh 



so|§:„-^| = 2"l|a„l|2/(r(rz)nr=iM- 

Next, set g{x) = [ILi ki|"'-"'exp (- ^"^, ai,„|x,|) /| J/(x)| in (l3iUD . Then 



/ n l^^r'"' I exp - XI (^i,n\Xi\ j dx 



vi=l / \ 1=1 



V/ ^'[ (f\\xi\^^^-Adn''-\x)dt 

■nil 2 Jo Ju\x\\„.„=t\ \7L-T / 



ll^"l|2 Jo J{\\x\\an=t} Vi = l 
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/"OO I' I ^ \ 

Jo illbiL =n IH- / 



|an||2 Jo ^{||^IU„=1} 

1 /• / " \ 9"ll/7 II r°° 

inlh y§"r^ \ :L, / r(n) n^i o^.n Jo 



Since 



lonih y§^-i \f^l * y r(n)nr=iai,nyo 

"^K" Vi=l / V i=l 

= 2" TT ( / f^'^'e-^^'^^dt J = 2" TT 



^i.n. I'^i I ^^ 



r(l + ai.„s) 






and 



/■OO 

/ t(^+^)"-^e~*rft = r ((1 + s)n) , 
Jo 
fl3XT]) follows. D 

Lemma 3.2. F2a; M > 1 anc? e G (0,1). T/ien i/iere exists a 6 > such that for all 
s G {—6, 6) \ {0} and every t G (0, M], we have 

(3.2.1) (1 - e)e-'^ < r(l + stY^' < (1 + e)e-*^. 

Proof. For s 7^ sufficiently close to 0, the term F (1 + s) can be estimated using L'Hopital's 
rule and the continuity of F'(a;)/F(x) at 1: 

(3.2.2) IhnF (1 + sf' = exp (^^^^ = e"^. 

Thus, we can choose 6 > such that for all s G {—6M, 6M) \ {0}, 

(3.2.3) (1 - £)^/^'e-T < F(l + sy/' < (1 + eY^^'e-^ 

Next, let s G (—5, 6) \ {0} and pick t G (0, M]. Using the change of variables s ^-H> st in (13.2.31) . 
we get (I3XT1) . D 

Lemma 3.3. Suppose there exists an M > 1 such that for all n E N, n > 2, and all 
i = 1, . . . ,n, ai^n < M. Then for every e > 0, there exists a 6 > such that for every 
s G [—6, 6) \ {0} and every n eN with n >2 we have 

(3.3.1) (1 - e)e-^ < H ^ (^ + "v"^) < (1 + ^)e"^- 

Proof. Fix n. Set t = ai^n in (13.2.11) . to get, for each i = 1, . . . ,n, 



;i_^)e-a.„7<r(l + sa,„)^/^<(l + 



e e 



-ai,n7 



Taking the product over i and using ^"^^^ aj_„ = n we obtain (13.3. ip . D 
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Lemma 3.4. For i = 1, . . . ,n let ti > 0. Subject to the restriction Yl^i=i ^i — "^^ ^^^ function 
f{ti, . . . ,t„) := nr=i ^i' achieves its global minimum when ti = t2 = ■ ■ ■ = tn = I, where it 
takes the value 1. 

Proof. Setting t** = 1 (by continuity) when ti = 0, the function /(ti,...,t„) := HILi^i' 
extends to the set where tj > for i = 1, . . . , n. Using Lagrange multiphers when all tj > 
and induction when one or several of the ti equal 0, the result follows. D 

Proof of Theorem \2.1[ Fix A; > and e G (0, 1) {k denotes a large parameter and e a small 
one). To obtain an upper bound for ri_„(x) = 'n-nr=i kil"*'" o^ ^a,7^' ^^ ^^e Chebyshev's 
inequality 



(3.4.1) PV\f{\^^r-'>t\<\E\i{ 
with s = s{k,e) > selected very close to 0, but fixed (so eventually s » n~^). Set 

(3.4.2) t = 2n^Elf[ 



It follows from fl3XT]) that 



n 



(3.4.3) ^:rMni^^r-<^^/^"|>i-^> 

so it suffices to prove 

(3.4.4) nt^/^o" < (1 + e)e-^ 

for some small sq (to be chosen below) and all n large enough; in particular, we always assume 
that s << 1, and that n » s"^ whenever both n and s appear in the same expression. Pick 
5 = 5{e) > such that (1 + (5)^ < 1 + e. By the choice of t (in 03.4.21) ) and Lemma [3?I1 we 
have 

Next we bound each of the three factors in the right hand side of ([3]) . From Lemmas 13.31 and 
13.41 we get, for all s > sufficiently small, 

(3.4.6) {\{^-^^i^Y\{l + S)e-'^. 

From the asymptotic expansion (I3.0.2p of F and the assumption n » s^^ we obtain 
(3.4.7) 

F in) \ "'- ^ /ex / (l + .)^/^- \ / ^ ^ / 1 \ \ "'-'' ^ ( 1 



F((l + s)n)y - \n) \{\^sY+^l^ ) \ \n ) ) - \n) \{\^sYI 
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Finally, by L'Hopital's rule, given any s > we have 
(3.4.8) lim (2n^)i/^" = 1. 

n—*oo 

Thus, we can select sq = so{e) > so small in f l3.4.7p that e/(l + sq)^^*" < 1 + 5. Choosing 
A^ = N{k, so{e)) such that for all n > A^ we have (2n^)^/'*o" <l + S, inequality (ISAij) follows. 
Observe that the hypothesis aj.„ < M on the renormalized weights entails that their geo- 
metric mean YYi=i 'h^n ^^ ^^^^ bounded above by M. To obtain 

Pr'{{^-er-VT<n\\\xr-\>l ^ 




2n^ 



basically all we need to do is to follow the same steps as before, but using s < instead of 
s > 0, and the preceding observation instead of Lemma [3. 4[ So we avoid the repetition. D 

The preceding proof works by respectively giving upper and lower bounds for nt^/^^ when 
s > and s < are close enough to zero. Since 

T [n) ^ 



nt^'''' = (2n'=)^/"" 



n 



T-r r (1 + a 






T{{l + s)n) 

and the first two factors on the right hand side approach 1 as ra -^ oo, concentration is 
controlled by the third factor. Furthermore, by Lemma [3.31 for all s sufficiently close to zero 
we have 

(3.4.9) (l-e)e-^< (nr(l + t 

so in order to determine how 




-|4 F (1 + aj^ns ) 

=1 ^i,". 



l/sn 



behaves it is enough to estimate its denominator. The concentration results in Theorems 12.51 
and 12.71 are obtained by giving sequences of weights for which the behavior of Y\a=i ^Tn i^ 
easily understood. 

Proof of Theorem \2.5[ Denote by j the integer part of n/{M + 1). Suppose first that 
n/{M + 1) is an integer. Then the sequence of weights aj^„ = M for i < j and aj.„ = 1/M 
for i > j satisfies ^"=0*^*." ~ ^- ^^ n/{M + 1) is not an integer, then jM + {n — j)/M < 
n while (j + 1)M + (n - j - 1)/M > n. Thus, there exists ate [1/M, M] such that 
redefining aj+i,„ = t (instead of 1/M) we have Yl^=o ^h"- ~ ^- Since for this family of weights 

T—rn a In M— 1 

lim„^oo nj=i '^i n" = ^ *^+^ 5 the result follows by using the same argument as in the proof 
of Theorem 12.11 D 

Theorem 12.71 is proven in the same way as the previous one, and Theorem 12. 8[ as Theorem 
12.11 save for the fact that we use the uniform probability P^~^ on the euclidean unit sphere. 
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and the corresponding expectation of YYi=i b^l^ o^^r §2 ^- This expectation is computed 
within the proof of Proposition 1 in [GluMij : 



(3.4.10) ^(ni^^i^)-X-(n''^'^)'''^" 



'rmV r(t) 



r(i) / r(i±^n)- 
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